(2) Japanese Patent Application Laid-Open No. 6-119308 (1994) 

"Simulation Method of Dynamic Characteristics of Semiconductor Device" 

The following is the extract relevant to the present invention: 

Fig. 1 is a connection diagram showing a simulation method of dynamic 
characteristics of a semiconductor device according to the preferred embodiment of 
the invention taking an insulated gate bipolar transistor as an example. As shown in 
Fig. 1, an insulated gate bipolar transistor (IGBT) 1 has a fundamental structure 
including a drain electrode 4 on the part formed by diffusing an n-type buffer layer 
and a p-type drain layer on a lower surface of an n-type substrate, a source electrode 
2 formed on the upper surface in a position to be contact with the source layer which 
is formed by diffusing a p-type channel formation layer and an n-type source layer 
and a gate electrode 3 formed on the source layer through a gate insulated film 5. 
Generally, the IGBT 1 is composite including a plurality of fundamental structures. 

A drain power source 1 1 (power source setting voltage Vd), a drain resistor 
10 and a drain inductance 9 are interposed between the drain electrode and the source 
electrode of the IGBT 1 to form an extemal circuit as a load circuit. A gate extemal 
circuit consisting of a gate power source 8 (power source setting voltage Vg), a gate 
inductance Ig as a wiring inductance 7 and a gate resistor rg as a wiring resistor 6 is 
interposed between the source electrode and the gate electrode. The IGBT 1 is 
generally switched on/off by the gate power source setting voltage Vd to control a 
drain current Jd as a load current. 

Next, assuming that short-circuit of the load is generated in the 
semiconductor device including the extemal circuit described above, a simulation 
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method will be described taking simulation of dynamic characteristics of a gate 
circuit including the short-circuit generated therein as an example. The simulation 
is performed by the combination of state equations (1) to (6) and an operation 



equation as follows: 
5 State Equations 

div(grad ^i^) = - (e/^) (Nd ~Na + p - n) ...(1) 

Dp = -(l/e)div JpH-G-U ... (2) 

Dn = + (l/e)div Jn + G-U ... (3) 

Jp = -eDp gradp - // p p grad y/ • • • (4) 

10 Jn = H-eDn gradn - // n n grad y/ ... (5) 

J = Jp + Jn - eT> grad y/ ... (6) 



Here, references represent the elements as follows: 
Operators 

D = d/dt : time differential operator 
15 div : scalar operator showing divergence 

grad : vector operator showing gradient 
Variables on calculation 

y/ : potential (voltage) J : current p : concentration of positive hole 
n : concentration of electron Jp : electron current 
20 Jn : positive hole current 

Constants to be provided during simulation 

e : charge of electron s : dielectric constant G : generation term 
U : recombination term Dp : diffusion constant of positive hole 
Dn : diffusion constant of electron /u p : mobility of positive hole 
25 // n : mobility of electron Nd : concentration of donor 
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Na : concentration of acceptor 
Operation Equation 

f {f (Vgm ~ Vg) dt + rg{q- go) - Ig (dq / dt) }dt = -lg(g- qo) 

Jq Jq 0 

Here, references represent the elements as follows: 
5 Ig : gate inductance rg : gate resistance 

qo : initial charge in gate insulated film q : charge as one parameter 
vg : gate power source setting voltage 
vgm : gate voltage as another parameter 

The wiring inductance Ig and the gate resistor rg of the gate external circuit 
10 are designated as circuit constants and the initial charge qo is designated per time 
variation. Further, the charge q as one parameter determined by the state equations 
is substituted into the operation equation, to thereby indirectly calculate the gate 
voltage Vgm as another parameter. 

Next, using a residual equation as follows as the modification of the 
15 operation equation described above, the charge q determined by the state equations is 
designated again and iterative calculations are performed according to Newton's 
iterative method until a residual f is fall within a predetermined range of convergence 
judgment error. 

/ (Vgm) = f {f (Vgm - Vg) dt + rg (q ~ qo) - Ig (dq / dt) J dt Ig (q - qo) 

Jq Jq 0 

20 That is, Vgm3 expected to have a value of "f (Vgm3) ^ 0 " can be 

obtained from the values of residuals f (Vgml) and f (Vgm2) to trail values Vgml 
and Vgm2 during Newton's iterative method. The residual f can fall within the 
predetermined range of convergence judgment error by repeatedly performing this 
step. 
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According to the method of preferred embodiment here, it is possible to 
check and adjust consistency between the state equations and the operation equation 
by the residual equation to implement calculation. Therefore, the amount of 
calculation errors within time variation can be considerably reduced. Further, as the 
5 calculation is shifted to the next time variation after the value of the residual f is 
reduced to fall within the predetermined range of convergence judgment error, it is 
possible to prevent vicious cycle of sequential accumulation of calculation errors 
extending through a large number of time variations accompanied by the progress in 
simulation and to improve accuracy in calculation. As a result, dynamic 
10 characteristics of the gate circuit including high-frequency vibration can be predicted 
with high precision. 

Further, as orders of both sides of the operation equation are equal, the value 
of the convergence judgment error s can be obtained from the following equation: 
s = f (Vgm) / lg(q - qo) 
15 It this value is assumed to 10*^ to lO""^, dynamic characteristics of the gate 

circuit including high-frequency vibration can be estimated with precision which is 
enough for actual use. 
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mt LX(D>y'- h«jEE5rvgm t i^tc t ^ m^m<Dmf^:^ 

[fell 



lg(dq /dt)o ) dt= -lg(q-qo) 



[If *^ 3 1 f (vgm) . ^^■rBm:^m^t l-c (12) 

^Srffll'X f (vgni)/lg(q-qo)= £ Xm^th^H^^^ii'^ 

mmA^ i =io-^~io-*o^iatciiR*s*-T?fi-^^Mt5 

t t 

f (Vgm) =J { J (V£in-VR)dt + rR('q-qo)- IfiCdq Alt)« } dt + lfifq-qo) 
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div (grad </> ) -= — (6/ t ) (Nd— Na+p 

Dp =-(l/e) div Jp +G-U 

Dn = + (l/e) div Jn +G — U 

J p = — e D p grad p — uP P grad </i 

J n = + e D n grad n — /in n grad (p 

J =Jp+Jn— fD grad </> 
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div : %iS;$r^-r^*7?t^^ 

grad : <fe)ia^^-#-^ii' hJl^i^M^ 40 

(p : :^T:y'>^/l^-t^£i:>hMBc. J : fl;^ifc^ P : iE7LjS 
n : J P : M^MM. J n : jEllMM 
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m^t^ -^mm (omwcx- 5 ^6 . ' > miu.mm v ^ u 

fzL V ^^r^ma ^ V(% Sr JlS^il ^6 b ^ - ; 

Vgm =Vg-rg(dq/dt) 
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qt^iK^:&S^(i) ~(6) ^M<^t\-2:^Mi^k.th^ 

^(7) :^mmt{:im'r^^mu^m^^/ux-\^^^(D 

[0 0 0 9] El 3 I G BT(D^nm^X-M^M 

fco fi5/^^SS»]^?F^Vgvfi. ^^^^tl5[H]£S<^^»'Y>' 
^^i^^^u-iyaz^jjmxnnwx^^. ^mi^mm 
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[0 0 0 6] MO S h Sr*t-^^3»^SB 

tLX<DmU'f-h^<^^^~'7h'7>'iy:^^ (IGB 
T) ^^'^O^'- b^ggi^y- h^>-^*>J^^>->;ig, 
<i:l>^^'- hSiargSr-^tf^-^. MHlgSco»j>m;^S^ 

vg fi>5r^- hm^E^^^iS^mJE. vgmti^^^^mi: l.X<D 
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t-. r<;:^^P^t::J:tb«. MO S h SrW-T^^^^g 

<Dm^mm(o^mm^i^'^ \.xmr^^m^m(oi&:^<Dmi\' 
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(^mi^jym^t'LxiiD^^m^^^^ :i tt-r^o 
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Ja{JoCVgiD-Vg)dt+rBCq Qo) -IgCdq /dt)o } dt= -lg(q-qo) 

[0 0 14] ^fefc. aMf(vgni) Sr^fa^:*ra^^ irM^m^ M-T^ t t^^o 

LT(l2)^Srffll/\ f (vgm) /lg(q-qo)= fX^^tL^ [00 15] 

iRm4^iJ^^^;i5 f =io-^~io-^o$5H(;iffi/h'r^^Xff-*40 [^6:4 ] 

f(Vsa)= fIr(Vgm-Vg)dt-hrg(q-qo)-lg{dq /dt)o } dt+lg(q-qo) 
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9 KV"i' W >'y^' v;^ 

10 K t-'-r i^m.^ 

11 Kw^'iiiig (fs^mmvd) 
Vgm <^'- hm& 
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